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ABSTRACT 

We revisit the kinematical data for 204 globular clusters in the halo of M87. Beyond 
3 r e gf along the major axis of the galaxy light, these globular clusters exhibit substantial 
rotation (~ 300±70 km s^ 1 ) that translates into an equally substantial spin (A ~ 0.18). 
The present appearance of M87 is most likely the product of a single major merger, 
since this event is best able to account for so sizable a spin. A rotation this large makes 
improbable any significant accretion of material after this merger, since that would have 
diluted the rotation signature. 

We see weak evidence for a difference between the kinematics of the metal-poor 
and metal-rich population, in the sense that the metal-poor globular clusters appear 
to dominate the rotation. If, as we suspect, the last major merger event of M87 was 
mainly dissipationless and did not trigger the formation of a large number of globular 
clusters, the kinematic difference between the two could reflect their orbital properties 
in the progenitor galaxies; these differences would be compatible with these progenitors 
having formed in dissipational mergers. However, to put strong kinematic constraints 
on the origin of the globular clusters themselves is difficult, given the complex history 
of the galaxy and its last dominant merger event. 

Subject headings: galaxies: individual (M87), galaxies: star clusters, galaxies: elliptical 
and lenticular, cD 
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1. Introduction 

Little is known about the kinematics of the 
halos of elliptical galaxies. Unlike spiral galaxies, 
ellipticals contain little gas which could serve as a 
kinematical tracer, and studies of the stellar light 
seldom extend further than 2 effective radii (r c g), 
due to the decreasing surface brightness and the 
difficulty obtaining spectra of the diffuse stellar 
light. Furthermore, companion galaxies, which 
could serve as kinematical probes in the outer 
halo, are usually not numerous enough and too 
distant to probe the range between 2 and 10 (or 
more) r c g. Recently, however, globular clusters 
and planetary nebulae have become popular test 
particles to probe the kinematics in the halo of 
giant ellipticals. With the commissioning of 8m- 
class telescopes, large samples of radial veloci- 
ties can be obtained, extending out to several r c g 
from the center of the galaxies. Globular clusters 
are especially numerous around central giant el- 
lipticals and therefore are probably the best trac- 
ers of the outer kinematics in these galaxies. 

The kinematics derived from the globular clus- 
ters can be used to constrain galaxy formation 
history. For example, cosmological N-body sim- 
ulations make specific predictions of the amount 
of spin that will result from the torques caused by 
the multiple accretion of objects (see Sect. 4.2); 
moreover, a formation via dissipational mergers 
predicts both varying kinematics of the globular 
clusters with radius and different kinematics be- 
tween the newly formed globular clusters and the 
ones brought in by the progenitors (see Sect. 4.3). 

Thus, the kinematics of a large sample of 
globular clusters provides not only constraints 
on the formation and evolution of the galaxy, 
but also a better understanding of the forma- 
tion of the globular cluster systems. In the cen- 
tral cD galaxy in Fornax (NGC 1399), the kine- 
matics of the globular clusters constrain the ori- 
gin of the large over-abundance of globular clus- 
ters (Kissler-Patig et al. 1998a). In the other 
dominant giant elliptical in Virgo (NGC 4472), 



the kinematics of the globular clusters suggest 
a formation via two massive gas-rich galaxies 
(Sharpies et al. 1998). 

Here, we revisit the central giant elliptical in 
the Virgo galaxy cluster, M87 (NGC 4486), and 
derive the kinematics in a region between 1 and 5 
r c g using published radial velocities for over 200 
globular clusters (Cohen & Ryzhov 1997). Pre- 
vious work on the kinematics of globular clusters 
around M87 also include studies from Huchra & 
Brodie (1987), Mould et al. (1990), and Brodie 
& Huchra (1991). These studies were based on 
20 to 45 objects. Interestingly, Mould et al. re- 
ported dynamically significant rotation in both 
M87 and M49. In M87 their study extends to 
~ 400" radius; they report a rotation along the 
major axis of 0.60 ± 0.27 km s _1 arcsec -1 , i.e., 
240 ± 108 km s _1 at their extreme radius. The 
extension of this result is the main motivation for 
this work. 

In the next section (Sect. 2) we briefly describe 
the sample and our new analysis before present- 
ing in Sect. 3 the velocity dispersion and rotation 
velocity both as a function of radius and an esti- 
mate for the dimensionless spin parameter A for 
the galaxy. We discuss our findings and present 
our conclusions on the formation history of M87 
in Sect. 4. 

2. Globular cluster sample and data anal- 
ysis 

2.1. The sample 

We used the sample of 230 globular clus- 
ter candidates around M87 for which Cohen & 
Ryzhov (1997) published radial velocities, includ- 
ing an additional cluster (ID 682) and the update 
on five velocities given in Cohen et al. (1998). 
The coordinates were obtained from Strom et 
al. (1981). Following Cohen & Ryzhov, we ex- 
cluded galactic stars by selecting a sub-sample 
of 204 globular clusters with velocities v r > 250 
km s _1 . We also excluded an additional candi- 
date at v r = 350 km s _1 . The velocity offsets 
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from the systemic velocity for the remaining ob- 
jects are well below the escape velocity of the 
galaxy. Also, we do not expect any torques or 
tidal influence on the velocities from the Virgo 
cluster, since M87 is located at its exact center 
(e.g., Schreier et al. 1982). The resulting sample 
is shown in Fig. 1, where we plot the position 
of all of the globular clusters. The symbol sizes 
reflect the difference between the cluster veloci- 
ties and our adopted systemic heliocentric veloc- 
ity for M87 of 1282 ± 9 km s _1 , taken from the 
RC3 (de Vaucouleurs et al. 1991). Open circles 
show approaching and filled circles show reced- 
ing globular clusters; the dashed line marks 1 
r cff (~ 100" , e.g., Goudfrooij et al. 1994) of the 
galaxy light and the solid line marks 3 r e ff . 

2.2. Analysis 

Cohen & Ryzhov (1997) comprehensively present 
and discuss the velocity dispersion of the globu- 
lar clusters as a function of distance from M87. 
The authors also discuss rotation, but only divide 
their data into two radial bins (within and out- 
side 180" ), and find no significant deviation in 
these two sub-samples compared with the rota- 
tion of the whole sample: about 100 km s _1 along 
a position angle (measured from north through 
east throughout this paper) of ~ 150°, defined as 
the angle of maximum positive rotation. 

We extend their study by investigating the ra- 
dial profiles of the rotation, position angle, and 
velocity dispersion. For each globular cluster, we 
use the nearest 75 data points in radius to mea- 
sure the above three parameters, allowing this 
sub-sample to drop to 30 data points at the ex- 
treme radii, 80" and 400". For each sub-sample, 
a maximum-likelihood fit determines the posi- 
tion angle, rotation amplitude, and dispersion at 
that radius (see Pryor & Meylan 1993, for exam- 
ple, for the basic approach). The fit measures 
the best sinusoid for the velocity data func- 
tion of position angle. We fix the mean veloc- 
ity of the sub-samples to be equal to the global 
mean velocity. The amplitude of the fit provides 



the rotation, the angle of maximum positive rota- 
tion provides the PA, and the standard deviation 
about the sinusoid provides the velocity disper- 
sion. The uncertainty in the fit is derived in two 
ways: using the classical covariance matrix and 
using a bootstrap. For the bootstrap uncertain- 
ties, at each location where a cluster exists, we 
draw a point from a Gaussian distribution with 
the mean given by the rotation amplitude and 
PA at that radius, and the standard deviation 
given by the dispersion at that radius. A velocity 
uncertainty of 50 km s _1 is also included. This 
sampling provides one realization and we simi- 
larly fit the three parameters. Generating 100 
realizations then allows for a distribution of val- 
ues at every radius; using the 16th and 84th per- 
centile values provide the 68% confidence band. 
The 68% confidence band and the classical la 
uncertainties are similar at each radius; in all of 
the figures we plot the confidence band based on 
the bootstrap technique. 

As mentioned above, each realization results 
from the assumption of a gaussian distribution in 
the velocity differences from the rotational veloc- 
ity. We have also used the actual differences for 
the boostrap simulations; i.e., we randomly draw 
with replacement from the distribution of differ- 
ences to generate a simulated dataset. There is 
no difference in the results using either a Gauss 
or the actual distribution. 

Cohen & Ryzhov (1997) quote a error of 100 
km s _1 for their velocities, but suggest that this 
is an overestimation; we therefore use an uncer- 
tainty of 50 km s _1 . However, we verirified that 
there is no difference in the results for either un- 
certainty, as is expected since the high velocity 
dispersion (about 370 km s _1 ) dominates the 
uncertainty in the estimation of the rotational 
velocity. 

Our results are shown in Fig. 2 and 3, where 
the projected rotation v • sin(i), projected veloc- 
ity dispersion a v , position angle, and v ■ sin(i) / a v 
(with their associated 68% confidence bands) are 
plotted against radius. Figure 2 shows the results 
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when the PA is allowed to vary and Fig. 3 shows 
the results for a fixed PA at f 20° (the mean value 
for the PA past 3 r e g). Figure 4 plots the individ- 
ual velocity measurements versus position angle 
in four radial bins of equal width (100") with 
mean radii of 1, 2, 3, and 4 r e ff, together with 
the derived rotation at these radii; rotation is 
clearly present in the outermost radial plot. Us- 
ing Monte-Carlo simulations and given a velocity 
dispersion of 370 km s" 1 , the probability of mea- 
suring a rotation this large, when no rotation is 
present, is less than 0.1%. 

3. Kinematics of M87 between 1 and 5 

r c flf 

Within 1 r e g, there is only marginal evidence 
for rotation in the sample, a result compatible 
with the lack of any significant rotation of the 
stars (Javis & Peletier 1991, Sembach & Tonry 
1996 and references therein). At ~ 1.5 r c ff, a 
group of globular clusters rotating along an axis 
offset by ~ 60° from the major axis possibly ex- 
ists; however, the result is only marginally signifi- 
cant. We will come back to this point in Sect. 3.2. 
The main feature is the monotonic increase of ro- 
tation past 2 r c fj. 

3.1. Kinematics beyond 2 r e ff 

Beyond 2 r c g, the globular clusters rotate 
roughly along the major axis of the galaxy (~ 
150°, e.g., Goudfrooij et al. 1994), with higher 
velocities in the SE. The position angle remains 
stable and the rotation amplitude increases from 
~ 50 km s _1 at 2 r c g to ~ 300 km s _1 at 4 
r c fj. The velocity dispersion remains constant at 
~ 370 km s~ x from 2 to 4 r c g. The increase 
in the velocity dispersion at the largest radii re- 
ported by Cohen &: Ryzhov (1997) is likely due 
to increasing rotational velocity. Accordingly, 
v • sin(i)/a v increases with radius from values 
around 0.2 to a value close to 0.8 at 4 r c g. Fur- 
ther evidence in support of the rotation is the 
correspondence of the major-axis position an- 



gle of the galaxy isophotes to the position an- 
gle of the maximum rotation. In addition, the 
amount of isophotal flattening is consistent with 
an isotropic oblate rotator given the amplitude 
of v/a. For an oblate rotator, we would expect 
an ellipticity of e ~ 0.25 for the v/a v ~ 0.5 seen 
at ~ 350", and e ~ 0.45 for the v/a v ~ 0.9 seen 
at ~ 400" (Binney & Tremaine 1987). 

McLaughlin et al. (1994) studied in detail the 
spatial structure of the globular clusters in M87 
and found the system to be elliptical, aligned 
along the major axis, with the ellipticity e in- 
creasing steadily to values of ~ 0.3 at about 300". 
The diffuse stellar light was studied to even larger 
radii (Liller 1960, King 1978, Carter & Dixon 
1978). These studies posit similar findings; the 
stellar light of M87 is flattened with an elliptic- 
ity smoothly increasing to values of e between 
0.35 and 0.4 at 350" - 400". These observations 
strongly support the rotation and its amplitude 
detected in the globular clusters at these large 
radii. 

In order to compare our rotation to various 
theoretical predictions in Sect. 4, we calculate 
the dimensionless spin parameter (Peebles 1971) 
A = J E 1 / 2 G' 1 M~ 5 / 2 , where J, E and M are 
the angular momentum, total energy, and total 
mass of the galaxy respectively, and G is the 
gravitational constant. A measures the total spin 
of a galaxy; however, since globular cluster radial 
velocities only exist over a limited radial range, 
we only measure a fraction of the total spin. We 
must therefore either extrapolate to include the 
whole galaxy or use the A' parameter introduced 
by Barnes (1992). A' is given by J / J maX) where J 
is the angular momentum of the subsample and 
>^max is the maximum angular momentum that 
the subsample can have. These two spin parame- 
ters are not directly comparable; i.e., a cold disk 
has A = 0.43 and A' = 1.0 (Barnes 1992). Be- 
low we calculate both parameters for comparison 
with theoretical predictions. 

To estimate the spin parameters, we depro- 
ject the velocity profile with various assumptions, 
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use the logarithmic potential given by Weil et 
al. (1998), and thus compute the mass, bind- 
ing energy, and angular momentum as functions 
of radius. For the logarithmic potential we use 
Vo = 550 km s -1 , R c = 100", and q = 0.85, cor- 
responding to a total mass of M87 of M(R = 
120kpc) ~ 1.2 x 1O 13 M . For the deprojection of 
the rotational velocity we use the Abel equation, 
and the surface brightness and density profile of 
Weil et al. Assuming that the galaxy is edge-on, 
we numerically integrate 

(1) 

to obtain the rotational velocity along the ma- 
jor axis. Where v is the density used in Weil et 
al., v$ is the internal rotational velocity, E is the 
surface brightness profile (given in Weil et al.), 
and v p is the projected rotation in Fig 3. The 
integral extends to r = oo, but our data only 
goes out to 400", thus requiring an extrapolation. 
We use four different extrapolations beyond our 
last data point: zero rotation outside of 400", ro- 
tation decreasing linearly with increasing radius 
(by 0.375 km s _1 arcsec ~ 1 ), constant rotation, 
and rotation increasing linearly with increasing 
radius (by 0.625 km s _1 arcsec _1 ). 

Two more assumptions are necessary to obtain 
the total angular momentum: the relation of the 
globular cluster's rotation profile to the galaxy 
rotation profile and the 2-D structure of the ve- 
locity field. We assume that the galaxy has the 
same rotation profile as the globular cluster sys- 
tem. For the velocity field, we assume constant 
rotation on cylinders and obtain the rotation at 
every position in the galaxy from the major axis 
rotation profile. The total angular momentum 
for the galaxy is then given by 

J(r) = 2tt [ r r' 2 dr' f v^{R) R p{R, z) sin6d9, 
Jo Jo 

(2) 



where R = r'sinO and z = r'cosO. Similarly, we 
substitute v max (R) for vAR) in equation 2 to ob- 
tain the maximum angular momentum, J ma x( r )- 
For v max (R), we use ^Jv^(R) + a 2 , where a = 

370 km s _1 , the measured projected velocity dis- 
persion. The calculated v max is very similar to 
the circular velocity at all radii, and using ei- 
ther one provides similar values for J max . J/Jmax 
equals A'. 

Figure 5 plots A' as a function of radius for the 
various assumptions. The confidence bands are 
calculated in the same manner as in Figure 4. 
In the extreme case, where the rotation veloc- 
ity drops to beyond our last observed point, A' 
rises steeply to values around 0.8 past 300". In 
the more realistic cases, where we assume that 
the points at the largest radii are near the maxi- 
mum net rotation and that the rotation curve re- 
mains approximately flat beyond our last point, 
A' reaches values of > 0.10 ±0.05. To decrease A' 
even further at r < 400", we would have to as- 
sume that the rotation velocity continues to rise 
steeply beyond 400" to values well above 700 km 
s _1 , an unrealistic probability. For any assump- 
tion, however, A' obtains a high value; putting all 
of the rotation at r < 400" results in A' around 
0.8, whereas allowing the largest rotation to be 
at r > 400" suggests a lower A' for the inner radii, 
but A' ~ 0.7 at larger radii. The significant ro- 
tation seen around r = 400" (v/a = 0.8) forces 
M87 to obtain high values of A' independent of 
the rotational-velocity extrapolation used. 

Equation 2 is used to calculate A' integrated 
from zero to R. However, A' can also be com- 
puted in radial bins (e.g. Hernquist & Bolte 
1992). Figure 6 plots the local A' integrated over 
a 40" region around R. 

Note that we make several important assump- 
tions. First, we speculate that M87 is seen edge- 
on, which could lead to an under-estimation of 
the derived angular momentum (by ~ 15% for 
i = 60°). However, the amount of flattening in 
the isophotes for M87 at large radii (e ~ 0.4) im- 
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plies that it is nearly edge-on, since any inclined 
orientation would require an unlikely more flat- 
tened system. Our second assumption is that the 
rotational velocity of the globular clusters repre- 
sents the rotation velocity of the total mass. If 
the globular clusters rotate twice as fast as the 
rest of the mass — an extreme case when com- 
pared to results for rotation velocities of galaxies 
in cosmological N-body simulations, see Sect. 4 
- we would overestimate J, i.e. A', by a fac- 
tor of two. Simulations, however, tend to show 
that the globular cluster system rotates more 
slowly than the stellar body (Hernquist & Bolte 
1994). Finally, we use a mass normalization for 
the logarithmic potential of Weil et al. (1998) of 
vq = 550 km s _1 ; using extreme values as low as 
vq = 400 km s _1 or as high as vq = 650 km s _1 
results in values of A' ~ 0.2 and ~ 0.05 respec- 
tively at 400". 

As well as A', we also estimate A. For A, 
we need the mass, angular momentum, and the 
binding energy. We cannot use A as measured 
from a subset of a bound system since A measures 
the total spin for a bound system; we must there- 
fore extrapolate to large radii to compare our 
measured A to in the theoretical predictions dis- 
cussed below. The theoretical calculations gen- 
erally use objects out to 1-2 half- mass radii to 
estimate the spin parameter. For M87, the half 
mass radii is anywhere from 150-300 kpc (assum- 
ing a distance of 16.5 Mpc), depending on where 
one truncates its halo. The globular cluster ve- 
locity data only extends to around 35 kpc, so we 
must make an enormous extrapolation to deter- 
mine the angular momentum. For the mass and 
binding energy, there is no extrapolation since 
the mass model incorporates the data at these 
radii. The integral over the density times the 
potential yields the binding energy, and the inte- 
gral of the density provides the mass profile. We 
extrapolate the net rotation by assuming that it 
is constant between our last point and the half- 
mass radius, and Equation 2 then provides the 
measurement of the angular momentum. With 



this extrapolation, A asymptotes to a value of 
0.18 at a radius around 600" (60 kpc) and stays 
constant beyond there. 

In summary, our best values for the dimen- 
sionless spin parameters of M87 are A' ~ 0.7 at 
r > 600" and A ~ 0.18 for the whole system. We 
will discuss the implications for the formation of 
M87 in Sect. 4. 

3.2. The kinematics around 1.5 r e g 

Two groups of clusters with 4-5 members each 
- immediately outside 1 r c g to the North and 
South of the galaxy center (see Fig. 1) — ap- 
pear to produce the rotation seen at ~ 1.5 r c g 
around an axis offset from the rotation axis at 
large radii. This pattern does not extend fur- 
ther out. The statistical significance is a little 
greater than 1 a and still compatible within the 
errors with a constant low rotation out to 2 r c g 
along the major axis. The offset axis may, how- 
ever, be due to incomplete azimuthal coverage of 
the data and, in particular, the lack of data in 
the NW around 2 r c gf. Alternatively, this rota- 
tion could be due to a separate group of clusters 
on orbit around M87 (e.g., from an accreted sys- 
tem). Abundances are available for 8 of these 
clusters (Cohen et al. 1998). The values scat- 
ter over 75% of the range spanned by the full 
sample (Mg b < 3A), neither supporting nor ex- 
cluding the accretion of a dwarf companion. The 
exact nature of these clusters remain to be inves- 
tigated, but do not influence our results beyond 
2 r eff . 

4. Constraints on the formation of M87 

We compare the kinematical properties of the 
globular clusters to galaxy formation scenarios, 
and briefly discuss why our results can not be 
explained by a single infalling satellite. From 
hierarchical formation scenarios the high value 
derived for the spin parameter suggests that the 
evolution of M87 was likely dominated by one 
major merger event. 
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4.1. Only an infalling satellite? 

We must verify that the measured rotation is 
characteristic of the system and not a group of 
globular clusters associated with a satellite or 
originating from a tidal tail of an interacting, 
stripped galaxy Evidence against these latter 
hypotheses is the alignment of the globular clus- 
ters with the major axis, the match with the stel- 
lar rotation, and the relatively smoothly increas- 
ing rotation and ellipticity of the whole system. 
Weil et al. (1998) recently reported diffuse (28 
mag arcsec~ 2 ) stellar light around M87 and sug- 
gest that it results from the accretion of a small 
spheroidal galaxy. The fan of light nearly aligns 
with the major axis with a large opening angle, 
suggesting that the orbit must have passed close 
to the center. However, their simulations show 
that the accreted galaxy must have been of low 
mass (10 9 — 1O 1O M0), and the associated number 
of globular clusters would be small. By compar- 
ison with numbers of globular clusters per unit 
mass for similar galaxies in the compilation of 
Zepf & Ashman (1993), the total number of glob- 
ular clusters associated with such a galaxy would 
be of the order 3 to 30; the number present in our 
magnitude limited sample would be a factor of 3 
to 4 lower. Thus, it is not possible for these ac- 
creted clusters to dominate the globular cluster 
system of M87 at 4 cS . 

Three other dwarf galaxies are potential can- 
didates for satellites: NGC 4486A, NGC 4486B, 
and IC 3443. They are projected within 10' of the 
center of M87 and within 30° of the PA of the ma- 
jor axis and could be at a similar distance as the 
central giant elliptical. Both NGC 4486A and 
NGC 4486B are counter-rotating with respect 
to the stars and globular clusters with velocities 
differences of ~ 800 and ~ 200 km s _1 , respec- 
tively, compared to the systemic velocity of M87. 
IC 3443 is a dwarf galaxy with Mb t — —15.5 
and it seems unlikely that it could have con- 
tributed many globular clusters to M87. Finally, 
we checked the abundances published by Cohen 
et al. (1998) for the clusters which dominate the 



rotation at large radii. No systematic pattern is 
visible; the values scatter over the range spanned 
by 75% of the range of the full sample, on the 
metal-poor side. 

In summary, no companion can significantly 
contribute to the globular cluster system of M87 
at ~ 4 c ff , so the rotation must be intrinsic to the 
globular cluster system associated with M87. 

4.2. An explanation for the high spin pa- 
rameter 

Peebles (1969, 1971) first estimated the result- 
ing angular momentum of a galaxy formed by 
gravitational instability. More recently (e.g. Ueda 
et al. 1994 and references therein), specific pre- 
dictions were made for hierarchical clustering 
models with large cosmological N-body simula- 
tions. The results of the different simulations 
agree well; all simulations appear to predict a 
lower dimensionless spin parameter A than we 
observe in M87. Values for A vary between 
0.01 and 0.1 with a mean around A = 0.05 and 
are mostly insensitive to cosmological parameters 
and fluctuation spectrum shape. In addition, A 
tends to be lower in high-density environments 
and appears to be anti-correlated with galaxy 
mass (Ueda et al. 1994, Efstathiou Sz Jones 1979, 
Barnes & Efstathiou 1987). 

The spin parameter that we derive for M87 is 
A ~ 0.18. The result is somewhat sensitive to our 
different assumptions, but should not be off by 
more than a factor of two (see Sect. 3.1). That is, 
the result for M87 is only marginally consistent 
with the simulations, and lies at the upper end 
of, or above, the simulation results. 

Interestingly, Warren et al. (1992) noticed 
in their simulations the apparent decrease of A 
at masses comparable to M87, which they at- 
tributed to their finite computational volume 
rather than to a physical process: the most mas- 
sive halos in simulations have no similar large 
halos to tidally torque them. That is, M87 
would have to have encountered a equally mas- 
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sive galaxy in order to gain enough angular mo- 
mentum to explain values of A around 0.1 or 
higher. High values of A have also been borne 
out in other types of simulations. Hernquist 
(1992, 1993) showed in merger simulations of 
two equal-sized galaxies how the spin parame- 
ter could increase in the outer halo. Hernquist k, 
Bolte (1992) made specific predictions for globu- 
lar clusters at different radii in the end-product 
of such a merger and noticed that A' rises from 
values around 0.05 in the center to values as high 
as 0.2-0.3 at several effective radii. 

These simulations seem to indicate that the 
most likely scenario is that M87 gained its large 
amount of angular momentum through a merger 
of two large galaxies of about equal mass. This 
scenario does not exclude a hierarchical forma- 
tion for the progenitor galaxies or further accre- 
tion of smaller companions (see Sect. 4.1), but the 
present appearance of M87 must have been cre- 
ated by one single major merger event. It seems 
very unlikely that several merging galaxies would 
fall into M87 on the same orbital path so as to 
not dilute each other's rotation signatures. Given 
the total stellar mass of M87 (several 1O 12 M ), 
the two progenitors must have had stellar masses 
around 1O 12 M or above, which effectively rules 
out that M87 is the product of the merger of two 
(several 1O 11 M ) spiral galaxies. The most likely 
progenitors are two giant ellipticals, perhaps cen- 
tral dominant galaxies in their own groups, given 
the large masses involved. In agreement with 
this, M87 shows no significant recent burst of star 
formation; the last major merger appears to have 
been gas-poor or mainly dissipationless, unless it 
happened well before z = 1. 

4.3. A comparison with NGC 4472 and 
NGC 1399 

We carried out a similar analysis for two other 
giant elliptical galaxies: NGC 4472 (M49), the 
brightest galaxy in the Virgo galaxy cluster, for 
which we used 57 globular cluster radial veloc- 
ities compiled by Sharpies et al. (1998); and 



NGC 1399, the central cD galaxy in the For- 
nax galaxy cluster, for which we used 74 globular 
cluster radial velocities compiled by Kissler-Patig 
et al. (1998). We only summarize the results here 
since we feel that the samples are too small to 
draw strong conclusions. However, we note that 
in both cases the rotation at large radii is signif- 
icantly smaller than in M87, resulting in smaller 
spins. 

In NGC 4472 the clusters out to ~ 1.5r cff have 
significant rotation, driving the value of A' up to 
~ 0.5, but it decreases to ~ 0.1 at the last data 
point (around 2.5 r e ff) and, extrapolating with 
constant rotation amplitude, A' remains approx- 
imately constant out to larger radii. The total 
spin parameter A around the virial radius is low, 
around 0.01. The result for A' is very sensitive to 
a particular group of ~ 10 objects at small radii. 

In NGC 1399 no significant rotation could be 
detected in the inner regions and only small ro- 
tation (~ 150 ± 100 km s _1 ) was measured past 
300" (~ 6r e ff). Accordingly, we derive values of 
A' around 0.3, but compatible with 0. The to- 
tal spin parameter A near the virial radius lies 
around 0.05. These derived values for the spin 
parameters are in good agreement with the val- 
ues obtained by cosmological N-body simula- 
tions for the formation of central cD galaxies. 

5. Differences between the red and blue 
globular clusters 

5.1. Kinematics of the globular cluster 
sub populations 

The globular cluster system of M87 is known 
from photometry to host two distinct population 
of globular clusters (Whitmore et al. 1995, El- 
son & Santiago 1996). In order to investigate 
whether these two sub-populations differ in their 
kinematics, we divide the kinematic sample ac- 
cording to the metallicity obtained spectroscop- 
ically for 150 candidates by Cohen et al. (1998). 
The "blue" and "red" populations are defined as 
globular clusters with [Fe/H]< —0.9 (80 candi- 
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dates) and [Fe/H]> —0.9 (70 candidates), respec- 
tively. This cut is motivated by the gap at similar 
metallicities in the globular cluster population of 
the Galaxy (e.g. Zinn 1985) and the metallicity 
corresponding to the gap in the V —I color distri- 
bution of the globular clusters in M87 (e.g. Whit- 
more et al. 1995, Elson & Santiago 1996). The 
analysis was performed in the same manner as for 
the full sample (see Sect. 2). Figure 7 shows the 
results using a fixed position angle (120°): veloc- 
ity dispersion a v and rotation velocity v-sin(i) of 
the blue (thick solid line) and red (thin solid line) 
samples are plotted function of radius. The 
68% confidence bands (dotted lines) are similar 
for both samples and are shown only for the blue 
sample. 

Blue and red samples have similar mean veloc- 
ities (1299 ±44 km s" 1 and 1292 ±47 km s" 1 re- 
spectively). The rotation in the red sample seems 
constant (and compatible with no rotation) at all 
radii, while the rotation of the blue cluster pop- 
ulation seems to increase with radius past 200" 
and is significant beyond 300". However, given 
the sample sizes, the rotation amplitudes of the 
two samples only differ by ~ la. 

The blue and red samples have similar veloc- 
ity dispersions at all radii, however, their rota- 
tion profiles differ, implying different total V 2 = 
(a 2 + v 2 ot ) profiles. Since both groups trace the 
same mass distribution, either their radial den- 
sity profiles are different (e.g., one population 
has a flatter density profile than the other) or 
their families of orbits are very different (e.g., 
the globular clusters of one population having 
either primarily tangentially or radially biased 
orbits). Strom et al. (1981) first noticed that 
the mean color of the globular cluster system in 
M87 becomes bluer with increasing distance from 
the center. Neilsen et al. (1998) suggest that 
the color variation is due to a changing number 
ratio of blue to red clusters with radius (simi- 
lar to what was seen in NGC 4472 by Geisler 
et al. 1996, and NGC 1380 by Kissler-Patig et 
al. 1997), suggesting different density profiles. 



At 400", the ratio of the total V 2 of the two 
samples is ~ 1.4. Assuming isotropic orbital dis- 
tributions, in order for the two populations to 
trace the same potential, this ratio translates into 
a ratio of ~ 1.4 for the exponents of the density 
profiles. This ratio does seem typical for the ex- 
ponents of the density profiles of the blue and 
red populations (e.g., the cases of NGC 4472 and 
NGC 1380) and for the ratio of the exponents for 
those globular cluster populations which follow 
the stellar light profiles and those whose density 
profile is more extended than the stellar light (see 
the compilation in Kissler-Patig 1997, for exam- 
pie). 

Similar velocity dispersion profiles for the red 
and blue populations are therefore compatible 
with the larger rotation of the blue globular clus- 
ters if the red population has a steeper density 
profile, which also explains the color gradient in 
the system. 

5.2. Constraints on the formation of the 
globular clusters 

M87 is known to host an extremely high num- 
ber of globular clusters with respect to its mass 
(e.g., McLaughlin et al. 1994 and references therein), 
as well as hosting at least two different globular 
cluster sub-populations. 

Several mechanisms can explain the presence 
of different globular clusters populations (e.g., 
Kissler-Patig et al. 1998b). In the case of M87, 
the last merger probably did not involve a large 
amount of star and globular cluster formation; 
a result supported by the findings of Cohen et 
al. (1998), who showed that the vast majority of 
globular clusters are old and must have formed 
in earlier events. In this case, the globular clus- 
ter kinematics cannot strongly discriminate be- 
tween different formation scenarios since they 
were dominated by this last merger. Indeed, 
there is only weak evidence that the blue globular 
clusters cause most of the net rotation. We also 
stress that the last major merger cannot explain 
the globular cluster over-abundance and point 
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to Harris et al. (1998) and Cote et al. (1998) for 
alternative scenarios in the case of M87. 

It is possible that the apparent concentration 
of the net rotation in the blue globular clusters 
may be a relic from the situation in the pro- 
genitor galaxies. Hernquist (1993 and references 
therein) showed that, in an equal-mass merger, 
the angular momentum of rotating components 
tends to be conserved. In the case of M87, the 
observations would be compatible with formation 
of the progenitors by dissipational mergers (see 
Ashman & Zepf 1992). The red globular clus- 
ters would have formed from the infalling gas 
that accumulated in the center of the progenitors 
and would have little angular momentum (e.g., 
Barnes & Hernquist 1996). Whereas the older 
blue globular clusters would have been partly 
spun up by the mergers which formed the pro- 
genitors of M87. But the orbital mixing caused 
by the last major merger event in M87 makes it 
difficult to draw firm conclusions on the orbital 
history of its globular clusters from their present- 
day kinematics. 

6. Summary 

We have re-analyzed published radial veloc- 
ities for 204 globular clusters around M87 and 
found significant rotation in the outer regions 
(> 250"). We derive a dimensionless spin pa- 
rameter for M87 of A ~ 0.18 (A' ~ 0.7) from 
the rotation of the globular clusters. A com- 
parison with cosmological N-body simulations, 
argues that such a high spin parameter is most 
likely the product of an equal-mass merger. A 
single major merger must have been the domi- 
nant event shaping the kinematics and appear- 
ance of the galaxy. There is some evidence that 
the rotation is confined to the metal-poor glob- 
ular clusters. If, as assumed, the last merger 
was mainly dissipationless, this kinematic differ- 
ence could reflect the situation in the progenitor 
galaxies of M87. 
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Table 1 

Rotation velocity and velocity dispersion as a function of radius 



Radius 


Rotation velocity 


Velocity dispersion 


Spin 


arcsec 


v ■ sin(i) [km s _1 ] 


a v [km s _1 ] 


A' 


83 


56 -38 


000+34 


n nc+0 06 
0.06_ ()04 


HO 


1 q+43 


1 Q+37 

al8_ 2 g 


n A/i +0.05 
°- 04 -0.03 


103 


1 Q+44 
18-18 


319±|| 


-°4I .03 


114 


-7-7+39 


326I 31 


n n/i+0-07 

°-04I .o3 


123 


/10+41 
42 -35 


ooc+29 
336l 28 


n nc+007 


1 O A 

134 


q-i +43 
8i -37 


ocn+30 
350I 29 


n no+0.06 


1 A 

14o 


71+53 


O/iQ+32 
oy:y -28 


U - U9 -0.06 


104 


cc+45 


qqc+30 
-JOU_ 26 


n iQ+0.07 
u - lz -0.06 


104 


RO +49 


/I £+25 

345_ 32 


n i/i+O.OS 


1 -7/f 

174 


CQ+SO 


Q r 7+33 

00 ' -30 


r, 1 /i +0.08 
°- 14 -0.08 


183 


qe+53 
JO -35 


c + 3 1 

362_ 32 


|-| 1 q + 0.08 

u - lo -0.08 


i no 
193 


/II +41 
41-41 


373_ 28 


|-| 1 q + 0.08 

u - lo -0.06 


ono 
zU3 


/II +45 
41-41 


369lg| 


n 1 n+0.09 
U - 1Z -0.06 


213 


/1C+41 

48I 48 


0-7C+26 

375I 34 


n 1 1 +0.10 
O-H-0.06 


223 




380+^ 


O-H+006 


233 




381±|| 


10+ 009 

u - iu -0.05 


244 


90 + 53 

^°-23 


375 + ^ 


09+ 09 

u - uy -0.06 


253 




366+£ 


09+ 09 

u - uy -0.07 


264 


451S 


3541^ 


07+ ' 10 

u - u '-0.07 


273 


451S 


qcQ+33 




283 


74 +51 
' 4 -68 


359t|l 


0-05lg; 05 


294 




362I27 


0-05l°; 05 


304 


134+53 


356+24 


07+ ' 12 

u - u '-0.07 


314 




3441- 


10+° 12 
u - lu -o.io 


324 


187±g| 


354±*» 


0-12±g;iS 


333 




355+2? 


o.i4tg:15 


344 


199±H 


362lf 3 


0.18l°;l! 


353 


99^+53 
^ z ' -63 


3681^ 


0.20l°;l« 


363 


218t«j 


3641^ 


20 +0 15 

u - zu -0.16 


373 


229±|1 


373+|? 


24 +0 ' 14 

u - zy: -0.15 


383 


267±£? 


381+tt 


31+0 12 

u - ol -0.12 


388 


3021?;? 


380l« 


34+0- 12 

u -°^-0.12 



All values were computed for a fixed position angle of 120° East 
of North. 

The radius is given in arcsec, e.g. for a distance to M87 of 16.5 Mpc, 
1 arcsec would correspond to ~ 80 pc. 

The velocity dispersion is corrected for rotation and velocity errors. 
The listed spin was computed in a 40" bin around each point 
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Fig. 1. — Position of the globular clusters around M87. The sizes of the symbols reflect the difference 
between the globular cluster velocity and the systemic velocity of the galaxy. Open circles represent 
approaching, filled circles represent receding globular clusters; the dashed and solid lines mark 1 r e g- and 
3 r e ff of the galaxy respectively. 
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Fig. 2. — Projected velocity dispersion, projected rotational velocity, position angle and v ■ sin(i)/a v as 
functions of radius. Dotted lines mark the 68% confidence bands. 
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Fig. 3. — Projected velocity dispersion, projected rotational velocity, and v • sin(i)/a v as functions of 
radius for a fixed position angle of 120°. Dotted lines mark the 68% confidence bands. 
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Fig. 4. — Velocities of the globular clusters plotted against their position angles for four radial bins 
centered on 1, 2, 3, and 4 r e g (assuming 1 r e g = 100"). The solid lines show the best fitted rotation. 
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Fig. 5. — The dimensionless spin parameter A' (integrated over radius) plotted against radius. The solid 
lines show A' for four different assumptions. In the extreme (unrealistic) case where we assumed the 
rotation velocity to drop to km s _1 beyond our last data point, A' increase rapidly to very high values 
beyond 200" . The three other cases — rotation velocity stays constant, increases slowly, or decreases 
slowly beyond our last data point — produce very similar values for A'. The dotted lines show the 68% 
confidence band in the case of a flat rotation curve beyond our last data point. The last data point 
of the sample is indicated as a vertical dotted line. The long dashed line shows the behavior of A' as 
extrapolated from our model. 
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Fig. 6. — Dimensionless spin parameter A' computed in 40" bins plotted against radius by assuming a 
constant rotational velocity beyond our last data point. 



19 



200 



100 





200 300 

radius (arcsec) 

400 | 1 1 1 1 1 — i 1 r 




200 300 

radius (arcsec) 

Fig. 7. — Projected velocity dispersion (upper panel) and projected rotation velocity (lower panel) plotted 
against radius for the blue ([Fe/H]< —0.9, thick lines) and red ([Fe/H]> —0.9, thin lines) globular clusters. 
The dotted lines show the 68% confidence bands for the blue population; the uncertainties for the blue 
and red samples are similar. 
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